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Abstract 

Both BaBar and Belle have found evidence for a non-zero width difference in the D°-D° 
system. Although there is no direct experimental evidence for CP-violation in D mixing 
(yet), we show that the measured values of the width difference y ~ Ar already imply 
constraints on the CP-odd phase in D mixing, which, if significantly different from zero, 
would be an unambiguous signal of new physics. 
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The highlight of this year's Moriond conference on electroweak interactions and unified 
theories arguably was the announcement by BaBar and Belle of experimental evidence 
for D°-D° mixing [TJ [2} E], which was quickly followed by a number of theoretical anal- 
yses H El El [71 El E]. While Refs. pQ El E] focused on the constraints posed, by the 
experimental results, on various new-physics models, Ref. [5] presented a first analysis 
of the implications of these results for the fundamental parameters describing D mixing. 
The purpose of this letter is to show that the present experimental results already imply 
constraints on a sizeable CP-odd phase in D mixing, which could only be due to new 
physics (NP). 

To start with, let us shortly review the theoretical formalism of D mixing and the 
experimental results, see Refs. [Till [TT] for more detailed reviews. In complete analogy to 
B mixing, D mixing in the SM is due to box diagrams with internal quarks and W bosons. 
In contrast to B, though, the internal quarks are down-type. Also in contrast to B mixing, 
the GIM mechanism is much more effective, as the contribution of the heaviest down-type 
quark, the b, comes with a relative enhancement factor (m 2 — m 2 d ) / (m 2 — m 2 ,) , but also 
a large CKM-suppression factor |V^ fe ^| 2 /|V^ s V^.* | 2 ~ A 8 , which renders its contribution 
to D mixing ~ 1% and hence negligible. As a consequence, D mixing is very sensitive to 
the potential intervention of NP. On the other hand, it is also rather difficult to calculate 
the SM "background" to D mixing, as the loop-diagrams are dominated by s and d 
quarks and hence sensitive to the intervention of resonances and non-perturbative QCD. 
The quasi-decoupling of the 3rd quark generation also implies that CP violation in D 
mixing is extremely small in the SM, and hence any observation of CP violation will be 
an unambiguous signal of new physics, independently of hadronic uncertainties. 

The theoretical parameters describing D mixing can be defined in complete analogy to 
those for B mixing: the time evolution of the D° system is described by the Schrodinger 
equation 

with Hermitian matrices M and T. The off-diagonal elements of these matrices, Myi 
and Ti2, describe, respectively, the dispersive and absorptive parts of D mixing. The 
flavour-eigenstates D° = (cu), D° = (uc) are related to the mass-eigenstates Di 2 by 

\D 1>2 )=p\D°)±q\D ) (2) 



with 



n\ M* — IF* 
H \ IV1 12 2 1 12 



(3) 



j>) M 12 -|r 12 ' 

|p| 2 + I q | 2 = 1 by definition. 

The basic observables in D mixing are the mass and lifetime difference of Dip, which 
are usually normalised to the average lifetime T = (Ti + T2)/2: 

_ AM _ M 2 - Mi _ Ar _ T 2 - Ti 
x = — - p , V=2f- ~^r- ■ (4) 
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In this letter we follow the sign convention of Ref. [5], according to which x is positive 
by definition. The sign of y then has to be determined from experiment. In addition, if 
there is CP-violation in the D system, one also has 

^ 1, <f> = arg(M 12 /r 12 ) ^ 0. (5) 

While previously only bounds on x and y were known, both BaBar and Belle have 
now found evidence for non-vanishing mixing in the D system. BaBar has obtained this 
evidence from the measurement of the doubly Cabibbo-suppressed decay D° — ► K + tt~ 
(and its CP conjugate), yielding 

y = (0.97 ±0.44(stat) ±0.31 (syst)) x 10~ 2 , 

x' 2 = (-0.022 ± 0.030(stat) ± 0.021(syst)) x HT 2 , (6) 

while Belle obtains 

y CP = (1.31 ± 0.32(stat) ± 0.25(syst)) x 10~ 2 (7) 

from D° —>■ K + K~ , tc + tt~ and 

x = (0.80±0.29(stat)±0.17(syst))xHT 2 , y = (0.33±0.24(stat)±0.15(syst)) x 10~ 2 (8) 

from a Dalitz-plot analysis of D° — > Kg7T + 7T~. Here yep y in the limit of no CP 
violation in D mixing, while the primed quantities x',y' are related to x,y by a rotation 
by a strong phase 6k*' 

y' = cos 5k-k — x sin 5kk, x' — x cos 5k* + y sin (9) 

Limited experimental information on this phase has been obtainted at CLEO-c [12] : 

cos 6 K * = 1-09 ±0.66, (10) 

which can be translated into 6k% = (0 ± 65)°. An analysis with a larger data-set is 
underway at CLEO-c, with an expected uncertainty of A cos 6k* ~ 0.1 in the next couple 
of years [13J; BES-III is expected to reach A cos 6k* ~ 0.04 after 4 years of running 
[14"] . The experimental result (fTUjl agrees with theoretical expectations, 6k-k = in the 
SU(3)-limit and 15^1^15° from a calculation of the amplitudes in QCD factorisation 
[T5] . Based on these experimental results, a preliminary HFAG-average was presented at 
the 2007 CERN workshop "Flavour in the Era of the LHC" p3]: 

x = (8.5^; 2 ) x KT 3 , y = (7.1± 2 ;°) x HT 3 . (11) 

Adding errors in quadrature, this implies 

- = 1.2 ±0.6. (12) 

y 
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The exact relations between AM, Ar, M 12 and T 12 are given by 



(AM) 2 - - (Ar) 2 = 4|M 



12 



ir 
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(AM)(Ar) = 4Re(M* 2 r 12 ) = 4|M 12 ||r 12 | cos< 



(13) 



Eq. (Ti"3"j) implies x/y > for |0| < 7r/2 and x/y < for tt/2 < \<f>\ < 3ir/2. In view of the 
above experimental results, we assume \<f>\ < 7r/2 from now on. 

As for the CP- violating observables, \q/p\ ^ 1 characterises CP- violation in mixing 
and can be measured for instance in flavour-specific decays D° f, where D° — > f is 
possible only via mixing. The prime example is semileptonic decays with 
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r(D c 



-X) - T{D° i+X) \q/p\ 2 - \p/q\- 



T{D° -> E-X) + T{D° i+X) \q/p\ 2 + \p/q\< 



(14) 



Although the B factories may have some sensitivity to this asymmetry, its measurement 
is severely impaired by the fact that D mixing proceeds only very slowly, resulting in a 
large suppression factor of the mixed vs. the unmixed rate: 



r(L» c 



x 2 + y 2 



T(D° -> i+X) 2 + x 2 + y 2 
Both in the K and the B system the quantity 



6 x 10 -5 . 



(15) 



A M 
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is very small, which however need not necessarily be the case for D's. From 
the general expression 



4 + r 2 + 4r sin . 
4 + r 2 — 4r sin . 



1/2 



(16) 
one derives 

(17) 



with r = |ri 2 /Mi 2 | and the weak phase defined in (JHJ). In the B system, one has r < 1 
(the current up-to-date numbers are r rs 7 x 10 -3 for Bd and r pa 5 x 10~ 3 for B s [16J), 
so that upon expansion in r 



1 + 



B,, 



12 



M 



12 



sin + 0(r^ 



Note that this formula refers to the definition <ft = arg(M 12 /r 12 ), which differs by +tt 
from the one used in Ref. [16J, <fi = arg(— Mi 2 /ri 2 ). For the K system, one finds r pa 
|Ar/AM| pa 2 from experiment, but now the phase 4> turns out to be small, so that 



1 + 



4r 



K 



4 + r' 



■0 + O(0 2 



1 + 



(19) 



In both cases, \q/p\ ~ 1 to a very good approximation. In the D system, however, there 
is no natural hierarchy r C 1, and of course one hopes that NP-effects induce \<f>\ ^> 0. In 
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Figure 1: \q/p\ 2 , Eq. ([20]) . as a function of the CP-odd phase 4> for the central experimental 
value f = 7.1/8.5. Solid line: full expression, dashed line: first order expansion around (ft = 0. 



this case, and because x and y have been measured, while \M\2\ and | Ti2 1 are difficult to 
calculate, it is convenient to express \q/p\o m terms of x, y, <fi, using the exact relations 
(115)1 . From (J3]) , and defining r = y/x, we then obtain 



V2{1 



f2) 



2(l + P) 2 + 16f^tan 



+8f tan (f)sec <p\ (1 + r 2 ) 2 — (1 — r 2 ) 2 sin (j) 



1/2 



(20) 



Note that for finite xy and = ±7r/2, diverges because — ► for <fi — > ±7r/2 from 
( TT3)) . In Fig. [1] we plot \q/p\ 2 as function of 0, for the central experimental value from 
HFAG, f = 7.1/8.5, Eq. (fTT)) . It is obvious that even for moderate values of <fi the small-0 
expansion is not really reliable. 

What is the currently available experimental information on CP- violating in D mixing, 
i.e. \q/p\ and <pl As already mentioned, the semileptonic CP-asymmetry (fill) has not 
been measured yet. What has been measured, though, is the effect of CP- violation on the 
time-dependent rates of D° — > K + tt~ and D° — > K~7r + . The BaBar collaboration has 
parametrised these rates as 



r(D°(t) -> oc e" 



•rt 



. x' 2 -I- v' 2 

R D + y/R^y'-Tt + - y - (rt) 2 
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r(D°(t) -> iT-7T + ) oc e~ r * 

and fit the _D° and D° samples separately. They find [2] 

y' = (9.8 ± 6.4(stat) ± 4.5(syst)) x 10~ 3 



(21) 



y'_ = (9.6 ± 6.1(stat) ± 4.3(syst)) x 10 



-3 



(22) 



Adding errors in quadrature, this means y' + /y'_ = 1.0 ± 1.1. BaBar also obtains values for 
x± which we do not quote here, because the sensitivity to the quadratic term in (1211) is 
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Figure 2: Left: y' + /yL as function of <j> for x/y = 1.2 (solid line) and x/y = {0.6,1.8} (dashed 
lines), from Eq. (fTT|) . Sk-k = 0. Right: y' + /y'_ as function of <p for x/y = 1.2 for 5^ = (solid 
line) and 8kk = ±65° (dashed lines). 



1 /2 

less than that to the linear term in y' ± . is the ratio of the doubly Cabbibo-suppressed 
to the Cabibbo-favoured amplitude, R X J? = \A(D° -> K + n-)/A(D° -> #-7T + )|. <W is 
the relative strong phase in the Cabibbo-favoured and suppressed amplitudes: 



A(D C 



K~ 



7T 



A(D( 



K+ir 



(23) 



the minus-sign comes from the relative sign between the CKM matrix elements V c a and 
V us . In the limit of no CP-violation in the decay amplitude, one has |v4(.D — > K~tt + )\ = 
\A(D — > K + 7i~)\, which is expected to be a very good approximation, in view of the fact 
that the decay is solely due to a tree-level amplitude. Then the relation of y'± to x, y and 
is given by 



{(y cos Sktt ~ x si n &Kn) cos 4> + (x cos 5k-k + V sin &k-k) sin 0} , 
{(y cos Sktt ~ x sin 5kk) cos 4> ~ ( x cos ^Ktt + V sin 5ktt) sin 0} . 



(24) 



Presently, the experimental result for y' + /y'_ is compatible with 1, although with con- 
siderable uncertainties. Any significant deviation from 1 would be a sign for new physics. 
In Fig. [2] we plot y' + /y'_ as function of cj), for different values of x/y and Skk- The figures 
clearly show that the value of y' + /y'__ is very sensitive to the phase 0, at least if Sk-k is not 
too close to —65°, which corresponds to the nearly constant dashed line in Fig. [2b. The 
reason for this dependence on Sxn becomes clearer if y' + /y'_ is expanded to first order in 



x(x 2 + 2y 2 ) cos 5k-k + y 3 sin Skit 



+ O(0 2 



(25) 



V- (x 2 + y 2 )(x sin 5 K ir - y cos 5k* 

For the central values of x and y, Eq. ( fTTI) . this amounts to 1 + 3.40 for 5k-k = 0, 1 — 3.30 
for 5ktt = +65° and 1 + 0.450 for 5k-k = —65°, which explains the shape of the curves in 
Fig. [2b. Evidently it is important to reduce the uncertainty of 5 Km which, as mentioned 
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Figure 3: Plot of |AT/Ar ; 



Eq. ([26]) , as a function of x/ysM and 0. 



earlier, will be achieved within the next few years. On the other hand, as shown in 
Fig. [2^, y' + /y'__, which depends only on the ratio x/y, but not x and y separately, is not 
very sensitive to the precise value of that ratio, but very much so to 0. The conclusion 
is that, even if x/y itself cannot be determined very precisely, y' + /y'_ will nonetheless be 
a powerful tool to constrain 0, at least once 8kw will be known more precisely. Already 
now very large values ~ tt/2 are excluded. 

Another, more theory-dependent constraint on can be derived from the value of y. 
This argument centers around the fact that (a) the experimental result (lllj) is at the top 
end of theoretical predictions ysM ~ 1% P2] and (b) new physics indicated by a non-zero 
value of always reduces the lifetime difference, independently of the value of x. This 
observation is similar to what was found, some time ago, for the B s system [18] . In order 
to derive it, we assume that new physics does not affect T^jl] so that Ti 2 = rf^. We then 
have 2|r 12 | = Ar SM and hence \ysu\ — IT^I/T. Using the relations (fT3"|) . we can then 
express the ratio |Ar/Ar | in terms of ysu, x an d 0: 



value of x (and any new physics in the mass difference). In particular one has y = 
for = ±7r/2 and i / 0, which follows from the 2nd relation f[T3"l) . Eq. (j2"6"|) is the 
manifestation of the fact that one does not need to observe CP-violation in order to 
constrain it. A famous example for this is the unitarity triangle in B physics, whose 
sides are determined from CP-conserving quantities only, but nonetheless allow a precise 
measurement of the size of CP-violation in the SM, via the angles and the area of the 
triangle. In Fig. (3] we plot |Ar/Ar SM | as a function of r = x/ysu- The zero at = ±7r/2 
is clearly visible. The experimental value \y/ysu\ — 0(1) then excludes phases close 
to ±7r/2. In order to make more quantitative statements, apparently a more precise 
calculation of y$M is needed. 

1 See, however, Ref. [19] for a discussion of the effect of tiny NP admixtures to Ti 2 . 




(26) 



This implies that new physics always reduces the lifetime difference, independently of the 




Two more CP-sensitive observables related to D° 
the Belle collaboration [3]: 
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K + K have been measured by 



x sin ( 



(27) 
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1 

2 



y cos ^ + 
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The present experimental value of yep is given in (JTj), that for Ar is (0.01 ± 0.30(stat) ± 
0.15(syst)) x 10~ 2 . Again, we can study the dependence of these observables on 0. In 
Fig. @Ji we plot the ratio ycp/y, which is a function of x/y and 0, in dependence on 0. As 
it turns out, this quantity is far less sensitive to <fi than y' + /y'_, the reason being that its 
deviation from 1 is only a second-order effect in 0: 



yep = y { 1 + 



4 i 2 2 4 

2 x + x y — y 



2(x 2 + y 



2\2 



+ O(0 4 ) 



(29) 



Hence, unless the experimental accuracy is dramatically increased, and because the results 
on y' + /y'_ and y/ysM already exclude a large CP-odd phase ~ ±tt/2, it is safe to interpret 
yep as measurement of y. In Fig. HJd we plot the quantity Ar/ycp- Also here there is a 
distinctive dependence on <p, with A^/y oc <fi for small 0, but the effect is less dramatic 
than that in y' + /y'_. 

In conclusion, we find that the experimental results on D mixing reported by BaBar 
and Belle already exclude extreme values of the CP-odd phase <fi close to ±7r/2. This 
follows from the result for y, which is close to the top end of theoretical predictions and 
can only be reduced by new physics, and from y' + /y'_ ~ 1. While y' + /y'_ — 1 vanishes in 
the limit of no CP-violation, y ~ Ar is a CP-conserving observable, which demonstrates 
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the usefulness of such quantities in constraining CP-odd phases. Also yep, A r and the 
ratio Ay/i/cp can be useful in constraining 0. As long as there is no major breakthrough 
in theoretical predictions for D mixing, which are held back by the fact that the D meson 
is at the same time too heavy and too light for current theoretical tools to get a proper 
grip on the problem, the long-distance SM contributions to x will completely obscure any 
NP contributions and their detection. The observation of CP violation, however, presents 
a theoretically clean way for NP to manifest itself and it is to be hoped that in the near 
future, i.e. at the B factories or the LHC, at least one of the plentiful opportunities for 
NP to show up in CP violation [20J will be realised. 
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